Scaling and Multiscaling Behavior of the Perimeter of Diffusion— Limited Aggregation 
(DLA) Generated by the Hastings— Levitov Method 
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In this paper, we analyze the scaling behavior of Diffusion Limited Aggregation (DLA) simulated 
by Hastings-Levitov method. We obtain the fractal dimension of the clusters by direct analysis 
of the geometrical patterns in a good agreement with one obtained from analytical approach. We 
compute the two-point density correlation function and we show that in the large-size limit, it agrees 
with the obtained fractal dimension. These support the statistical agreement between the patterns 
and DLA clusters. We also investigate the scaling properties of various length scales and their 
fluctuations, related to the boundary of cluster. We find that all of the length scales do not have 
a simple scaling with same correction to scaling exponent. The fractal dimension of the perimeter 
is obtained equal to that of the cluster. The growth exponent is computed from the evolution of 
the interface width equal to j3 = 0.557(2). We also show that the perimeter of DLA cluster has an 
asymptotic multiscaling behavior. 

PACS numbers: 64.60.al, 05.20.-y, 61.43. Hv, 68.35.Fx 



I. INTRODUCTION 

Diffusion- limited aggregation (DLA), introduced by 
Witten and Sander [1] , has been shown to describe many 
pattern forming processes including dielectric breakdown 
[2], electrochemical deposition [3, 4], viscous fingering 
and Laplacian flow [5] etc. 

This model begins with fixing a seed particle at the cen- 
ter of coordinates in d dimensions. By releasing random 
walkers from infinity and allowing them to stick as soon 
as they touch the cluster, a fractal pattern grows. 
This procedure is equivalent to solving Laplace's equation 
outside the aggregated cluster with appropriate bound- 
ary conditions. The walker sticks to a point on the sur- 
face of the aggregate with a probability proportional to 
the local field strength at that point (the harmonic mea- 
sure). 

In two dimensions, since analytic functions automat- 
ically obey Laplace's equation, the theory of conformal 
mappings provides another mechanism for producing the 
shapes. This method has been directly used by Hast- 
ings and Levitov (HL) to study DLA [6]. These authors 
showed that DLA in two dimensions can be grown by 
using successive iterating stochastic conformal maps. In 
the present paper, we are interested in these off-lattice 
DLA patterns generated by this method. 

We present some evidence that the patterns generated 
by HL method have the same statistics as DLA clus- 
ters simulated according to the original definition. In 
the first part of the paper, we calculate the fractal di- 
mension of the cluster patterns by direct measurements. 
We use two different methods, first, using the scaling rela- 



tion between the average gyration radius of the generated 
patterns with their size, and the second, calculating the 
density two-point correlation function. We show that the 
results agree with the fractal dimension of DLA clusters. 

In the second part of the paper, we investigate the 
scaling properties of various length scales and their fluc- 
tuations, related to the boundary of the patterns. We ex- 
amine whether they follow a simple scaling relation with 
a same correction to scaling exponent, or their scaling 
behavior is governed by the multiscaling property. 

The multiscaling of DLA clusters, proposed by 
Coniglio and Zannetti [7], stands for space dependent 
fractal dimension according which a whole set of scaling 
exponents exists. It has been also claimed by Somfai, 
et.al. [8, 9, 10], that these scaling claims are misled by 
finite size transients, and DLA obeys simple scaling and 
all length scales scale with the same fractal dimension. 

However our simulation for clusters generated by HL 
method, shows that the growth exponent defined by the 
interface width, differs from the fractal dimension, and 
we find no correction to scaling exponent for it. Fur- 
thermore, we extend the concept of multiscaling to the 
boundary of the clusters and we find that the asymptotic 
behavior of the boundary also agrees with the multiscal- 
ing property. 



II. THE HASTINGS LEVITOV METHOD 

In the quasi-stationary approximation, the probability 
density of finding a particle satisfies the Laplace equation 
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where the zero boundary condition on the boundary of 
cluster 9C, implies the sticking of the particle upon ar- 
rival, and the later condition states that ijj(z) is indepen- 
dent of any direction at infinity. 

The probability of cluster growth at a certain point z of 
the boundary of the cluster is determined by the har- 



monic measure 



dP(z) = \Vif>(z)\dl, 



(3) 



where dl is a boundary element containing the point z. 

According to the Riemann mapping theorem, there ex- 
ists a conformal map that maps the exterior of the unit 
circle to the exterior of the cluster. Hastings and Levi- 
tov constructed this map using the iteration of conformal 
mapping [6]. The function <fi\^(w) maps the unit circle 
to a circle with a bump of linear size \f\ at the point 
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FIG. 1: A circle in the w-plane is mapped to a C„ in the 
z-plane by $ n (w). The same function maps a circle with a 
bump at # n +i to a C n +i- 



^e(w) = e l9 ^ (e- w w). (5) 

The parameter < a < 1 determines the shape of the 
bump, for higher a the bump becomes elongated in the 
normal direction to dC, e.g. it is a line segment for a = 1. 
In this paper we set a = g for which the bump has a 
semi-circle shape. 

A cluster C n consisting of n bumps can be obtained by 
using the following map on a unit circle 

$„(» = <px u8l ° <P\ 2 ,e 2 ° • • • ° (j) Xn .0 n (w), (6) 

which corresponds to the following recursive relation for 
a cluster C n +i (see Fig. 1), 

$„+iH = $n(<l>\ n+1 ,e n+1 (w)). (7) 

Since z = (fr n (w), one can obtain that 

dl = \<p' n (e* e )\d9, (8) 

where the prime denotes for differentiation. 
In order to have fixed-size bumps on the boundary of the 
cluster, since the linear dimension at point w is propor- 
tional to | <& Tl (""j) | 1 , one obtains 

Ari+1 = K(i+or (9) 

From Eq. 2 and Eq. 8 can be obtained that 

dP = \Vi)\\<5>'\d9 = d6, (10) 

indicating that the numbers 9 n have a uniform distribu- 
tion in the interval 0^9^ 2ir. 



In this paper our analysis is based on the boundary 
of the clusters and we need to have a uniform data on 
the boundary. This can be done formally by using a 
uniform series of {/? s } s=1 during the conformal mapping 
from a unit circle to the boundary of the cluster i.e., 

{w s = e 1,3s } s . This procedure can not be applied op- 
erationally, because in order to have a reasonable data in 
the fjords, one has to set S ^> n which needs very long 
simulation time. 

Barra et al, [11] have focused on the branch points of the 
map and introduced another approach for selecting the 
series {{3 S }- Following their approach, we define and 
as " Right" and " Left" branch points of the function 
4>\ n fi n in the following map, respectively 

where \a„ — a^|/27r is the fraction of the unit circle cov- 
ered by the bump. Each new bump creates two new 
branch points on the boundary and in case of probable 
overlapping with previous branch point, some of the older 
ones will be removed. So the maximum number of branch 
points will be 2n. If w k ' be a branch point of the fcth 
bump without overlapping by the next (n — k) bumps, 
it would be an exposed branch point of the map but 

the pre-image of the branch on the unit circle will change 
r r.l , R,L 

from w k to w k n 

L ) = $„(<; L ), (12) 

such that 

w k,n = 4>lt,e n ° ' ' ' ° 4>\l hl ,B h+1 ( w k' L ) ■ (I 3 ) 
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The solvability of Eq. 13 determines whether the branch 
point remains exposed, and then by mapping them one 
gets a reasonable image of the fjords. 



III. SIMULATION 

The simulation of the boundary of DLA clusters of 
different sizes is carried out using the algorithm discussed 
in the previous section. We set the parameter a = \, for 
which the function 4>\,0( w ) is analytically invcrtiblc. 

At the nth step, 9 n and A n are determined as follows. 
9 n is selected from a uniform distribution in the range 
[0, 2tt], and then A„ is computed using the Eq. 9. Af- 
ter determination of As and 9s and computing exposed 
branch points Wu'^i together with Eq. 6, the boundary 
of each cluster is determined. 

We generated 2000 clusters of number of bumps 10 3 ^ 
TV ^ 5 x 10 4 and 200 clusters of TV = 10 5 . A typical 
growth cluster is shown in Fig. 2. All average quantities 
which will be discussed later are taken over the simulated 
cluster ensemble. 



IV. DIRECT CLUSTER ANALYSIS 

In this section we do some direct measurements based 
on the geometry of clusters obtained from simulation. 
These include computation of the fractal dimension of 
generated DLA clusters and size-dependence of the vari- 




FIG. 2: Boundary of a typical simulated DLA cluster consist- 
ing of TV = 10 s bumps generated by using the HL algorithm, 
with a = | . The plotted shell is used to study the muliscaling 
properties of the boundary in Sec. V. The width of the shell 
is magnified by a factor of 10. Inset: a close-up view of the 
cluster. 



ance of gyration radius of the clusters. We find a good 
agreement between our results and ones obtained from 
the analytical approach in [12, 13]. We also measure the 
density correlation function— which, to our knowledge, 
has not been computed yet for the HL method— and we 
investigate its dependence on the size of the cluster. We 
find that the large-size behavior of the function corre- 
sponds to an expected correlation exponent a which is in 
a good agreement with the computed fractal dimension. 



A. Scaling of Gyration Radius for DLA Cluster 

The fractal dimension D c of DLA clusters generated 
by HL method has been previously computed from the 
Laurent expansion of the conformal map, cf. Eq. 6, equal 
to D c = 1.713(3) [12, 13]. The error in the last digit is 
indicated in parentheses. This has been obtained from 
the scaling relation between the first coefficient of the 
Laurent series of (j> n {w) and the size of DLA cluster. 
Since the first coefficient is proportional to the radius 
of the cluster, this motivates us to measure the fractal 
dimension directly using the scaling relation between the 
average gyration radius R g of the cluster and the number 
of bumps —or cquivalcntly the cluster size— TV, i.e., R c g ~ 
N Uc : where v c = 1/D C . 

The result is shown in Fig. 3(a). We find that v c = 
0.581(2), in good agreement with previous results. 

Another important result pointed out in [14] is the 
sharpness of the distribution of the first laurent coeffi- 
cient. It has been shown numerically that the rescalcd 
distribution width of squared first-laurent coefficient 
tends to zero as TV goes to infinity. Here, we check the 
same idea for the gyration radius of the clusters. The 
standard deviation of gyration radius is calculated from 

(T c = (Rg 2 ) — (Rg) 2 , where (•) denotes the ensemble 
average over simulated clusters of size TV. The rescaled 
<7 C as a function of TV is plotted in Fig. 3(b). As can 
be seen from this figure the fluctuation tends to zero for 
larger cluster size. This suggests that the rescaled distri- 
bution function of gyration radius of the clusters tends 
asymptotically to a 8 function. 

In order to investigate the asymptotic scaling behavior 
of er c , we proceed in the same way as [8, 10], where the 
authors suggest that all of the length scales I in DLA 
have a scaling relation with TV, like 

£~ N 1/D (a + bN-' fi ) (14) 

with a single universal exponent tp = 0.33(6). 
Our computation shown in Fig. 3(c) agrees with this scal- 
ing relation but with a different exponent of (p = 0.45(5), 
indicating that in the limit TV — * oo, the fluctuation of 
gyration radius has an asymptotic scaling behavior as 
that of the gyration radius, and nevertheless, the expo- 
nent seems not to be universal (this will be confirmed 
again in the following section for other length scales). 
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FIG. 3: (a) The average gyration radius of clusters Rg, versus 
the number of bumps N. (b) Rescaled standard deviation of 
gyration radius a c / R c g versus N. (c) Rescaled standard devi- 
ation of gyration radius (a c /N 1/D , D = 1.711) versus N~ v . 
The error bars are almost in the same size as the symbols. 

B. Density Correlation Function 

In this subsection, wc compute the two-point correla- 
tion function c(r), defined as 

where p(r) is density at position r, and the average is 
taken over all the points that belong to the cluster. For 
isotropic clusters the density correlation depends only on 
distance r. 

For self-similar fractals, c(r) should have the scaling form 
of c(r) ~ r~ a , where the exponent a is named co- 
dimensionality and is equal to a = d— D c , where d is the 
embedding dimension. 

Operationally, we proceed as follows to determine the 
function c(r) . For each sample in the ensemble of clusters 
of a fixed size, we cover the cluster by a two dimensional 
square lattice. Then for each lattice site belonging to the 
cluster, we consider an annulus around it with mean ra- 
dius of r and thickness of a lattice spacing. The density 
of the cluster points in the annulus is then proportional 
to the two-point correlation function at distance r. The 
average is then taken over both all lattice points in the 
cluster and all clusters in the ensemble. This procedure 
is repeated for annulus of different mean radius. 
We find that, for intermediate distances, the function 
c(r) exhibits a power-law behavior with an exponent a 
depending on the cluster size N. This behavior is shown 
in Fig. 4 for three different sizes. The values of the ex- 
ponent a as a function of the inverse size of the cluster is 
depicted in the inset of Fig. 4. In order to determine the 



value of the exponent in the large-size limit, we fit a poly- 
nomial curve to the data. We find that it extrapolates to 
a = 0.29(1), whose value is checked not to be affected by 
the degree of the fitted polynomial. This value is in good 
agreement with the aforementioned relation a = d — D c , 
with D c ~ 1.71. 
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FIG. 4: Two-point density correlation function c(r) for three 
different cluster sizes N. The difference in the slope of the 
solid lines indicates the size dependence of the correlation 
exponent a. The graph for N — 10 5 is shifted downward by 
0.12. Inset: the exponent a versus 1/N. The solid line is a 
polynomial fit of order 5, which yields the asymptotic value 
of q = 0.29(1). The error bars are almost in the same size as 
the symbols. 



V. BOUNDARY ANALYSIS 

In this section, we study the scaling properties of var- 
ious length scales related to the boundary of DLA clus- 
ters produced by HL method. We find that the fractal 
dimension of the boundary is the same as the DLA clus- 
ter, in agreement with the same conclusion reported in 
[15], for DLA clusters produced according to the original 
definition. We also check the simple scaling relation Eq. 
14 for various length scales including the gyration radius 
R b g , maximum radius R ma x and width R w of the bound- 
ary and their fluctuations. We find that all these length 
scales do not obey the scaling form Eq. 14 with a single 
exponent tp. 

Finally, wc check the multiscaling hypothesis for the 
boundary of the clusters and we will present evidence 
pointing to the existence of such anomalous scaling. 
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A. Scaling of boundary characteristic lengths 

Each cluster boundary is divided into segments such 
that ith segment has a length l i: and the distance of 
the midpoint of the segment from the center of mass is 
denoted by Ri. During the calculations, this procedure 
attributes a weight of U to each distance Ri and measures 
the following length scales in a more delicate manner. 



1. Gyration radius of the boundary, R[ 



The gyration radius of the boundary R b g is defined by 

Rg = j- hR 2 , where L is the total length of the 

boundary, and the sum runs over all segments on it. The 
fractal dimension of the boundary Db can be measured 
by using the scaling relation R b g ~ L Vb , where Vb = 1/Db- 
Fig. 5 shows the ensemble average of gyration radius 
versus the average length of the boundary. We find that 
Vb = 0.587(4). This indicates that within the statistical 
errors, a DLA cluster generated by HL method and its 
boundary have a same fractal dimension i.e., v c = Vb- 
This result is the same as one obtained before for DLA 
patterns grown according to the original definition [15]. 
It may be considered as another evidence that the pat- 
terns generated by method of iterated conformal maps 
proposed by Hastings and Levitov agree statistically with 
ones originally introduced by Witten and Sanders. We 
have also checked the scaling of R g with the cluster size 
and we found the same behavior as R g with N. 
The inset of Fig. 5 shows the plot of the rescaled stan- 
dard deviation of R b i.e., ^/N 1 / " against N~ v . We 
find that tp = 0.31(5), in agreement with Eq. 14. 



2. Maximum radius of the boundary, Rmax 

The other length scales we discuss here, are the lengths 
related to the maximum value of Ri in each cluster 
boundary represented by Rmax in Fig. 6. We observe 
from Fig. 6(a) that the ensemble average of R m ax scales 
with size N, with v max = 0.571(1), different from the 
gyration radius exponent. As shown in Fig. 6(b), the 
rescaled R max follows the simple scaling behavior of Eq. 
14, with a quite different exponent of tp = 0.18(5) from 
the proposed universal value of tp = 0.33(6) in [8, 10]. We 
also checked this simple scaling behavior for the rescaled 
standard deviation of R ma x, in agreement with Eq. 14 
(see Fig. 6(c)). 
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FIG. 5: Average gyration radius of the boundary R^, versus 
the average length of the boundary L. Inset: rescaled stan- 
dard deviation of the gyration radius {o b / N 1 ^ D , D = 1.711) 
versus . The error bars are almost in the same size as 
the symbols. 



T Si h(Ri — R) 2 , where the mean radius of the cluster 
is R= iE,W.- 

The growth exponent (3 can be obtained from the evo- 
lution of the interface width R w ~ N 13 . As shown in 
Fig. 7(a), we obtain the growth exponent for DLA clus- 
ters generated by HL method equal to j3 = 0.557(2). 
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3. Interface width, R w 

According to the analogy between the DLA grow- 
ing cluster and non-Euclidean growing interfaces, 
the interface width R w can be defined by R w = 



FIG. 6: (a) The ensemble average of the most-far boundary 
segment from the seed Rmax, versus the number of bumps 
N. (b) Rescaled R max (i.e., i? ma;c /iV 1/D , D = 1.711), ver- 
sus 7V~^. (c) Rescaled standard deviation of R ma x (i.e., 
a max /N 1/D , D = 1.711) versus N~ v . The error bars are 
almost in the same size as the symbols. 
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We checked the correction to scaling for the exponent, 
according to Eq. 14 represented in Fig. 7(b), and we 
conclude that no correction exists. The fluctuation of 
the interface width (see Fig. 7(c)) exhibits a simple scal- 
ing relation of form Eq. 14, with a correction to scaling 
exponent of tp = 0.58(5). This exponent is very different 
from those obtained for hitherto mentioned length scales, 
and far from its proposed universal value. 

The scaling properties of the interface width, appar- 
ently deviates from the simple scaling of Eq. 14, which 
has been proposed in [10] on refuting the multiscaling 
property of DLA cluster. The deviations of these bound- 
ary related length scales from the simple scaling behavior, 
motivated us to check an extension of the multiscaling 
property (previously applied for the mass of DLA clus- 
ters) to the length of the perimeter of clusters. 
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FIG. 7: (a) The average interface width of the boundary of 
cluster Rw, versus the number of bumps N. (b) Rescaled 
interface width (i.e., R w /N 1/D , D = 1.711) versus N~ v . 
(c) Rescaled standard deviation of the interface width (i.e., 
a w /N 1/D , D = 1.711) versus N~ v . The error bars are al- 
most in the same size as the symbols. 



B. Multiscaling analysis of the boundary of DLA 
clusters 



In this section, we extend the concept of multiscaling, 
previously used for the mass of the DLA clusters [16, 
17, 18, 19], to the length of the border of DLA. Our 
measurement for the perimeter of DLA clusters of size 
up to 10 5 particles (or bumps), reveals the multiscaling 
behavior of the border. 

For each cluster size, we generated an ensemble of DLA 
clusters by using the HL method and the perimeter of 
each sample has been determined as described in Sec. II. 
We proceed as follows: for each sample perimeter in the 
ensemble of size N and average gyration radius of R b g , a 



FIG. 8: Examples of the density profile of the boundary length 
of the clusters within a shell of rescaled radius x = r/R b g , 
represented for four different sizes. Inset: log-log plot of the 
density profile at a certain rescaled radius of x = r/R b g = 
0.8 (the dashed line in the main figure). The slope of the 
fitted solid line yields D(x = 0.8) = 1.67 for N = 100000. 
This figure summarizes the procedure we applied to obtain 
the functions D(x) in Fig. 9. The error bars are almost in 
the same size as the symbols. 



shell of radius r and of width dr (which is about the linear 
size of a bump) is drawn (see Fig. 2 for illustration). 
Then we measure the density profile g (r, R g ) defined as 



g (r, R g ) dr = dl, 



(16) 



where dl is the total length of the boundary within the 
shell of radius r. 

The plot of g (r, R g ) as a function of the rescaled ra- 
dius x — r/Rg, within 0.1 < x < 2, is shown in Fig. 
8 for four different sizes. This function has a maximum 
for distances around the gyration radius of the cluster. 
Assuming the scale invariance of the density profile [17], 
the multiscaling exponent D(x) can be defined as 



g(r,R g ) = C{x)R^-\ 



(17) 



where C{x) is a scaling function. Thus, the multiscaling 
exponent can be obtained using the following relation 



D(x) = 1 



d In g{r,R g ) 



d\nR„ 



(18) 



The inset of Fig. 8 shows the procedure we used to 
determine the multiscaling exponent as a function of x. 
At each x, the values of the density profile are read from 
Fig. 8 for each cluster size of gyration radius R g , and 
then D(x) is determined by Eq. 18. 

The whole behavior of D{x) for different size intervals 
is shown in Fig. 9. This shows that the function D{x) 
does not tend to a constant value as size increases, and 
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FIG. 9: Multiscaling fractal dimension D(x) of the boundary 
for different cluster sizes as a function of x = r /R h g . 



there is a maximum around x ~ 1.2 whose location does 
not depend on the size of cluster. Using the curves of 
Fig. 9 (and other similar curves obtained for other clus- 
ter sizes which not shown in the figure), we also estimated 
the value of D(x) at each x in the limit of N — > oo. As 
shown in Fig. 9, D(x) is not constant and varies with 
x, suggesting a multiscaling behavior. We therefore con- 



clude that the perimeter of the DLA clusters generated 
by HL method does not have simple scaling, and thus a 
set of scaling exponents is needed to be described. 



VI. CONCLUSION 

We studied scaling properties of DLA clusters gener- 
ated by the Hastings-Levitov method. First, we calcu- 
lated the fractal dimension of the clusters by direct ana- 
lyzing of the DLA patterns in agreement with the previ- 
ous results. We also computed the two-point correlation 
function of the mass of the cluster, and we found that 
in the large-size limit, it agrees with the obtained fractal 
dimension. 

In the second part of the paper, we focused on the bor- 
der of the DLA clusters and we investigated their scaling 
properties. We found that the fractal dimension of the 
perimeter is equal to that of the cluster. We checked the 
simple scaling behavior for various length scales including 
the gyration radius, maximum radius and the interface 
width of the boundary, together with their fluctuations. 
We found that all of these length scales do not have a 
simple scaling with a universal correction to scaling ex- 
ponent. The growth exponent has been obtained from 
the evolution of the interface width. Finally, we found 
that the perimeter of DLA displays an asymptotic mul- 
tiscaling property. 
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